Quantum phase transition in quantum wires controlled by an external gate 
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We consider electrons in a quantum wire interacting via a long-range Coulomb potential screened 
by a nearby gate. We focus on the quantum phase transition from a strictly one-dimensional to a 
quasi-one-dimensional electron liquid, that is controlled by the dimensionless parameter nxo, where 
n is the electron density and xq is the characteristic length of the transverse confining potential. 
If this transition occurs in the low-density limit, it can be understood as the deformation of the 
one-dimensional Wigner crystal to a zigzag arrangement of the electrons described by an Ising order 
parameter. The critical properties are governed by the charge degrees of freedom and the spin sector 
remains essentially decoupled. At large densities, on the other hand, the transition is triggered by the 
filling of a second one-dimensional subband of transverse quantization. Electrons at the bottom of 
the second subband interact strongly due to the diverging density of states and become impenetrable. 
We argue that this stabilizes the electron liquid as it suppresses pair-tunneling processes between 
the subbands that would otherwise lead to an instability. However, the impenetrable electrons 
in the second band are screened by the excitations of the first subband, so that the transition is 
identified as a Lifshitz transition of impenetrable polarons. We discuss the resulting phase diagram 
as a function of nxo. 

PACS numbers: 71.10.Pm, 64.70.Tg, 75.40.Cx 



I. INTRODUCTION 

Electron correlations in one-dimensional (ID) systems 
are especially pronounced due to the restricted avail- 
able phase space. Most prominently, the quasi-particle 
concept for the electron liquid, that accurately describes 
metals in higher dimensions, breaks down in ID giving 
rise to a correlated Luttinger liquid with the concomi- 
tant spin-charge separation. At strong interactions, the 
formation of Wigner crystals and commensurate Mott in- 
sulating states are expected. Technological advances in 
the fabrication of quantum wires and carbon nanotube 
systems with a high tunability nowadays allows the con- 
trolled study of such strong correlation phenomena^— 
see Ref. for a recent review. 

While transport and spectroscopic properties of quan- 
tum wires in the Luttinger liquid regime have been 
the subject of much attention over the years both 
experimentally 2 ^ 4 ^— and theoretically^^ - — much less is 
known outside this regime where a description solely in 
terms of a Luttinger liquid becomes insufficient. This 
is, in particular, the case close to quantum phase tran- 
sitions at which the number of Fermi points changes 
and, consequently, a Fermi energy Ep vanishes. In the 
approximation of a non-interacting electron gas, these 
transitions are directly reflected as rounded steps in the 
conductance ) 14 i 15 G(Ep,T), that sharpen up as the tem- 
perature T decreases. The influence of electron-electron 
interactions close to these transitions is, however, only 
partially understood*^— It is important to realize that 
the limits Ep — > and T — > in the approach to the 
quantum phase transition do not commute, and, in par- 
ticular, the critical conductance G(Q, T) is governed by 



the underlying quantum critical points As a Fermi en- 
ergy Ep vanishes at the transition, one necessarily enters 
a regime where temperature T is larger than character- 
istic energy scales of the electron liquid and, as a result, 
phenomena beyond a Luttinger liquid description become 
important. Close to the quantum phase transition, one 
can distinguish two characteristic energy scales, namely 
the Fermi energy Ep for charge excitations and the spin 
exchange constant J for spin excitations. The intermedi- 
ate temperature range Ep 3> T J, where the charge 
sector still sustains plasmon excitations while the spin 
sector is already incoherent, has caught some attention 
recentl y 26 ' 33 ' 34 It was pointed out by Matvee-^ that 
in this spin-incoherent regime the conductance should 
qualitatively deviate from the prediction for the non- 
interacting electron gas, which might explain additional 
features in the conductance of quantum wires observed 
experimentally 1 ' 18 

In the present work, we consider the quantum phase 
transition from a strictly one-dimensional electron liquid 
to a quasi-onc-dimensional state with the aim to iden- 
tify the nature, i.e., the universality class, of the tran- 
sition at T = 0. In the absence of any interaction, the 
electrons form subbands whose separation is controlled 
by the confining potential. The quantum wire is one- 
dimensional if only the lowest subband is occupied and 
undergoes a transition to a quasi-one-dimensional state 
by filling the second subband. Without interaction, this 
is just a Lifshitz transition where the number of Fermi 
points increases from two to four. It is the purpose of 
this work to investigate how this transition is modified 
in the presence of Coulomb interaction between electrons 
that is screened by a nearby gate. 

For spin-polarized electrons this question was already 
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FIG. 1: Transition from a one-dimensional to a quasi-one- 
dimensional state, (a) At weak coupling, nan S> 1, the tran- 
sition is triggered by filling a second subband upon tuning the 
chemical potential /i2 through zero, (b) At strong coupling, 
nan -C 1, the transition corresponds to the deformation of a 
ID Wigner crystal to a zigzag configuration. 

addressed in a series of works by some of the authors^ - — 
It was found that in the limit of weak interactions, 
na,B 3> 1 j where n is the one-dimensional electron density 
and as is the Bohr radius, the transition is still a Lifshitz 
transition but in terms of polarons, i.e., dressed fcrmionic 
quasi-particles. As the electrons start occupying the sec- 
ond subband, see Fig. [TJa), they propagate with a very 
small velocity. The relatively fast Luttinger liquid fluc- 
tuations within the filled first subband can follow these 
electrons adiabatically thus dressing them with a cloud 
of plasmons. It was shown by Balcnts^ that the residual 
couplings of the dressed fermions to the Luttinger liq- 
uid are irrelevant at the transition. A peculiarity of this 
transition concerns processes where pairs of electrons are 
transfered between the two subbands^S. Although such 
processes do not influence the nature of the transition, 
they determine in fact the ground state for a finite den- 
sity of polarons, i.e., pair-tunncling is dangerously irrel- 
evant at the critical point. As polarons start populating 
the second subband, the system can gain energy by pair- 
tunncling and, as a result, a BCS-likc gap opens up as 
a secondary effect. As a result of the dangerously irrele- 
vant pair-tunneling operator, the Lifshitz transition thus 
separates two CI phases in the notation of Ref. HH that 
possess only a single gaplcss charge mode. 

As the interaction strength increases, the universality 
class of the transition, however, changes. For strong in- 
teractions, na,B *C 1, a one-dimensional Wigner crystal 
forms, and the transition corresponds to the splitting of 
the crystal into two rows with a zigzag arrangement of 
the electrons, see Fig. [TJb). This transition can be de- 
scribed by an Ising order parameter. Correspondingly, 
the transition is now of Ising type and separates two CI 
phases with a single gapless phonon mode. The residual 
coupling between the Ising critical degrees of freedom 
and the plasmon excitations leads to logarithmic correc- 
tions to Ising criticality, that were analyzed in Ref. Hl|. 
Although this coupling turns out to be marginally irrel- 
evant, the renormalization group (RG) flow generates an 
enhanced SU(2) symmetry as the velocities of the plas- 
mons and the Ising critical degrees of freedom approach 
each other in the low-energy limit. However, at the same 
time the mean velocity decreases resulting, e.g., in a di- 
vergence of the specific heat coefficient as a function of T. 



The multicritical point separating the Lifshitz transition 
at weak and the Ising transition at strong coupling has 
not been identified so far. 

In this work, we extend the analysis of these previ- 
ous works to electrons that are not spin polarized and 
investigate the influence of the spin degree of freedom 
on these transitions. One-dimensional two-subband sys- 
tems of electrons have been investigated before by many 
authors i 35 " 41 For a related theoretical study in the con- 
text of cold atomic systems see Ref. Of particular in- 
terest in the present context are the works of Varma and 
Zawadowski^ that was motivated by the physics of fluc- 
tuating valence compounds, and of Balents and Fisher— 
that analyzes the phase diagram of the two-chain Hub- 
bard model. We discuss the relation of our results to 
these earlier works in detail in the main text. The zigzag 
transition of the Wigner crystal was also studied in the 
context of ion traps^r— as well as dipolar cold gases4£ 

Experimentally, the influence of the confining poten- 
tial (and, thus, the interaction strength) on the transi- 
tion from a one-dimensional to a quasi-one-dimensional 
state was investigated by Hew et al~ with the help of 
conductance measurements on a weakly-confined wire. 
The absence of the first quantized plateau in a range of 
confining potential strengths was interpreted as a signa- 
ture related to the formation of a two-row Wigner crystal 
state. 

The organization of this article is as follows. In section 
[IT] we specify the model used in our study and construct 
a mean-field phase diagram. The transition for strong 
interactions, na,B *C 1, is discussed in section HTT1 while 
the weak coupling regime, nas 3> 1, is addressed in sec- 
tion |IVJ We conclude in section IVl with a summary and 
discussion of our results. 



II. MODEL OF THE QUANTUM WIRE 

In the following, we specify the model of a quantum 
wire created by gating a two-dimensional electron gas 
(2DEG). The 2DEG forms at the interface (chosen to be 
in the xy-plane) between two different semiconductors 
which provide a steep confining potential in the perpen- 
dicular direction due to their band structure mismatch. 
Applying a bias voltage to a metallic split gate at a dis- 
tance d from the 2DEG leads to an additional confining 
potential V con f(y) that restricts the motion of the elec- 
trons in the y-direction, thus creating a wire along the 
x-direction. In addition, the gate screens the Coulomb 
interaction between the electrons. 

The system is described by the model Hamiltonian 

H = f+V coni + V int , (1) 

where the single-particle part comprises the kinetic en- 
ergy, T, and the confining potential, V con { , whereas Vint is 
the Coulomb interaction energy. We work in units h = 1 , 
Hb = 1, and 47reo = 1. 
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FIG. 2: Mean-field phase diagram of the Hamiltonian ([TJ) 
as a function of the inverse one-dimensional electron den- 
sity n _1 and the oscillator length xq — l/^ 



'mfi, both mea- 
sured in units of the Bohr radius as = e/(e 2 m). As xo in- 
creases, a transition occurs at the (red) solid line from a one- 
dimensional to a quasi-one-dimensional state. At low densi- 
ties l/(nas) 3> 1, this transition corresponds to the deforma- 
tion of a ID Wigner crystal into a zigzag configuration. At 
high densities, l/(nas) <S 1, the transition is triggered by the 
filling of a second subband. The thick dotted line indicates 
where the interaction energy equals the subband separation 
so that the band picture ceases to be well-defined. For large 
xo, the two-dimensional limit is approached. For a derivation 
of the various lines and regimes, see text. 



Assuming a parabolic confining potential characterized 
by the frequency f2, the single-particle terms read 



(2) 



where the index i labels the electrons and m is the effec- 
tive electron mass. The interaction between electrons is 
given by 



V, 



int 



(3) 



with the screened Coulomb potential 

e V|r| V 1 " 2 + (^) 2 



(4) 



where e is the electron charge and e is the dielectric con- 
stant of the host material. The mirror charge in the 
gate ensures that the potential falls off as a dipolc field 
U(r) w 2e 2 d 2 /(e\r\ 3 ) for large distances, r > d. 



1. Characteristic length scales 

The Hamiltonian ([lj is characterized by four character- 
istic length scales, (i) The mean particle spacing is given 
by the inverse one-dimensional density, 1/n, i.e., the 
number of particles in the quantum wire per unit length. 
(ii) The extension of the system in the lateral direction 
is quantified by the oscillator length, x$ = l/v^mfi. Fi- 
nally, the electron-electron interaction is characterized by 
(Hi) the Bohr radius, = e/(e 2 m), and (iv) the distance 
to the gate, d, specifying the strength and range of the 
interaction potential (Q}, respectively. In the following, 
we assume that the Bohr radius is much smaller than the 
distance to the gate, as -C d, as is usually the case for 
quantum wires; for example, for GaAs based systems one 
can estimate as — 10 nm and d > 100 nmi2 

To understand the interplay between the different 
length or energy scales, a mean-field phase diagram of the 
model as a function of the length scales 1/n, xq, as, and 
d can be constructed by minimizing the dominant terms 
in the Hamiltonian in various regimes and comparing the 
associated ground state energy scales. In particular, one 
can distinguish the cases where either the single-particle 
part of the Hamiltonian dominates over the interaction 
energy V int (Sec. Ill A 2|) or vice versa (Sec. Ill A3|) . By 
comparing the Fermi energy, Ep ~ n 2 /m, which is the 
characteristic energy for single-particle physics, with the 
typical interaction energy at distances of order of the 
mean particle spacing, U(l/n) ~ ne 2 /e, one finds that 
the two cases are distinguished by the ratio of the mean 
particle spacing and the Bohr radius. For high densities, 
nas > 1; the single-particle energy dominates whereas 
for low densities, naB < 1, the interaction energy domi- 
nates. 



2. Single-particle limit: multi-subband quantum wire 

In the high density limit, naB 3> lj interactions are 
weak. Thus, it is appropriate to first chose the eigen- 
basis of the single-particle part of the Hamiltonian con- 
sisting of T + V con f, see Eqs. @, and then analyze the 
effect of interactions in that basis. The single-particle 
cigenbasis is given by product-wavefunctions of traveling 
waves along the wire and the oscillator eigenfunctions in 
the transverse directions. The system is one-dimensional 
(ID) as long as only the lowest oscillator level is occu- 
pied. By comparing the Fermi energy with the oscillator 
frequency f2, one obtains the condition 



x < 



(5) 



A. Mean-field phase diagram 

Before analyzing the Hamiltonian ([TJ , let us discuss its 
qualitative features. 



that is shown as a solid line in Fig. [2] The transition 
from a single- to a two-subband quantum wire that is 
at the focus of this work occurs for nxo ~ 1. Upon 
further increasing the density n or relaxing the confining 
potential, i.e., increasing x , more and more subbands 
are populated. 
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The oscillator levels are well defined as long as the os- 
cillator frequency Q is larger than the typical Coulomb 
energy U(l/n), which translates to the condition xq < 
\ja,Bln, see dotted line in Fig. [2] For even larger oscil- 
lator length Xq one crosses over into a two-dimensional 
regime where the subbands are washed out. 



get filled and analyze the effect of interactions. While a 
single subband can be described as a Luttinger liquid dis- 
playing spin-charge separation, the interactions between 
the subbands couple spin and charge degrees of freedom. 
This leads to interesting modifications as compared to the 
spin-polarized case ; 29 ' 31 that will be discussed in Scc. lIVI 



3. Interaction dominated regime: Wigner crystal 

For small densities, nas *C f, interactions dominate 
over the kinetic energy of the particles, and the system 
behaves almost classically. In that case, the potential en- 
ergy, V^ on f + Vint, should be minimized first, resulting in 
a Wigner crystal state of electronsj 47 ' 48 As a function of 
oscillator length xq, there is a competition between the 
interaction and the confining potential. At strong con- 
finements, xq — > 0, a one-dimensional Wigner crystal is 
the ground state of the potential energy T4 on f+Vi n t. If the 
confinement is relaxed, a transition to a zigzag Wigner 
crystal occurs, followed by subsequent transitions to more 
two-dimensional multirow structures ! 29 ' 30 ' 49 " — Compar- 
ing V con { with the energy gain obtained by relaxing the 
Wigner crystal from a strictly one-dimensional arrange- 
ment, U(r = (l/n,0)) — U(r = (l/n,y)), one arrives 
at the stability criterion for the one-dimensional Wigner 
crystal state, 

indicated by the solid line in Fig. [5] Note that, at very 
low densities, n < 1/d, the screening of the interaction 
becomes important and modifies the transition line be- 
tween the one-dimensional and the zigzag crystal. At 
even lower densities, n < as/d 2 , the Wigner crystal melts 
(not shown in Fig. [2]). 



III. WIGNER CRYSTAL QUANTUM WIRE 

At low electron densities, nas <C I, the interaction 
energy dominates over the Fermi energy, see Fig. [2j and 
the classical ground state is a good starting point to de- 
scribe the physics of the system. Minimizing the po- 
tential energy V con f + Vim 01 Eq. (p} one arrives at a 
Wigner crystal state for the electrons^ - At low elec- 
tron densities or strong confinement such that Eq. (JB]) 
is obeyed, this Wigner crystal corresponds to a one- 
dimensional arrangements of electrons with equilibrium 
position (x'j,y l j) = (ja,0) of the j th electron, where 
a = n~ x is the lattice spacing. If the confinement is 
relaxed and the oscillator length xq increases, a transi- 
tion from a one-dimensional arrangement of electrons to 
a zigzag Wigner crystal takes place, see Fig. [IJb). In 
the following, we consider the properties of this quantum 
phase transition. 



A. Charge sector 

The excitations of the one-dimensional Wigner crystal 
correspond to fluctuations of the electrons around their 
equilibrium positions, 



B. Outline 

The focus of this work is the evolution of the quan- 
tum wire from the strictly one-dimensional limit to a 
quasi-one-dimensional state. At fixed Bohr radius as 
and distance to the gate d, a transition to a quasi-one- 
dimensional state occurs upon increasing oscillator length 
xq or density n, see Fig. [2] This transition is the sub- 
ject of the present paper with a focus on the interplay 
of charge and spin degrees of freedom. In the strongly 
interacting Wigner crystal regime, the increase of n or xq 
leads to a transition from a ID to a zigzag Wigner crystal. 
Charge and spin degrees are governed by very different 
energy scales so that, as a starti ng poi nt, they can be 
treated separately as done in Refs. I29ll3fl (charge) and [HI 
(spin). The question as to how the spin sector affects 
the Wigner crystal transition is addressed in Sec. IIII1 In 
the weakly interacting regime, the increase of n or xq 
results in the subsequent population of subbands. We fo- 
cus on the transition where the second subband starts to 



(u 



X ^-3, y A, (7) 
a a / 



which arc the phonon modes. Deep in the ID regime, 

the only low-energy mode is the Goldstone mode ujj cor- 
responding to a translation of the Wigner crystal along 
the wire direction, (u X j, u y j) = 0). As the confining 
potential is weakened and the quantum phase transition 
is approached, there appears another low-energy mode 
u^. It corresponds to the alternating oscillation of elec- 



trons in the direction transverse to the wire, (u x 



i) 



u^(0, (— 1) J ). The remaining two phonon modes, the lon- 
gitudinal out-of-phase phonon mode i4l((— 1) J 0) and the 
transversal phonon Uq (0,1), have gaps determined by the 
Coulomb repulsion and the confining potential, respec- 
tively. Close to the transition both these gaps are of the 
order Q. 

The Lagrangian density C = C\\Q+£± n +£i n t governing 
these low-energy excitations in the continuum limit is 
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given by 



C\\ = 



m 
~2v 

m 
2ri 



c mt = a a 



(a T ^) 2 + ^ (a^) 2 + r + s {u^Y 

(8b) 
(8c) 



X LtQ I , 



where r is the imaginary time. The longitudinal veloc- 
ity W|| is determined by the compressibility, mv 2 /n = 
d 2 E /dn 2 , where E is the ground state energy; in 



the limit d 



< n < a 



B , it evaluates to v 2 Q 



2n/(m 2 as) ln(8nd)j22 The control parameter r in (]8b[) 
changes sign as function of Xo, namely r = l/m 2 (xo~ 4 — 
x c~ 4 ) with x 0c = (2e/(7C(3)me 2 n 3 )) 1 / 4 . The transverse 
velocity vj_ n and coefficient of the quartic term s read 
v\t, = n/(m 2 a B ) In(2) and s = 93C(5)n 4 /(8™ 2 a B ), re- 
spectively. Finally, the coupling between the two modes 
is given by A = 21C(3)e 2 n/(4e). 

Away from the transition r > 0, the transverse mode 
uj^ corresponds to an optical phonon with minimal fre- 
quency oc y/r. With increasing oscillator strength xq, this 
gap decreases and vanishes at a critical value xq c , signal- 
ing the instability of the one-dimensional Wigner crystal. 
For r < the phonon field condenses with a finite ex- 
pectation value corresponding to a zigzag Wigner crystal 
phase, see Fig. [TJb). For r = 0, the Lagrangian C±_ v 
corresponds to a critical one-dimensional Ising model 
(s > 0). Neglecting the coupling between modes, A, the 
transition from the ID to the zigzag Wigner crystal is, 
thus, in the Ising universality class and the transverse 
mode has a gap which scales linearly with the distance 
from the transition point, A^J: = r. If the Wigner crystal 
was pinned and the positions of particles along the wire 
were fixed, the transition would break the reflection sym- 
metry in the confining plane. However, the presence of 
the mode wjj, i- e -j the fact that the crystal may deform in 
the longitudinal direction, makes the zigzag order non- 
local^ 

It is well-known that the transverse field Ising model 
can be alternatively represented in terms of a fermionic 
degree of freedom \1/ with Lagrangian density^ 



£ fcrm = ^g T y 



V_Lt, 

2 



($d x iS> + h.c.) + rtf** . (9) 



The longitudinal plasmon itjj couples to the most relevant 
operator (u^) 2 ~ vfitvj/ f the Ising model so that the 
8cl) can be rewritten in the fermionic 



interaction term 
formulation as 



ferm 



A 



(a^o) 



(10) 



where, for simplicity of notation, we suppressed in 
Eqs. © and (jTU)) rcnormalizations of coupling constants. 

The model £|| + £ { "™ + £ l ™ t m an d its critical proper- 
ties were analyzed and discussed in Ref. l3lL It was found 



that the critical renormalization group flow of the model 
parameters depends on the ratio of velocities, uj_?rA'l|o- 
If u^tt < ui|o, which is the case for quantum wires (see 
above) , the interaction A is marginally irrelevant and de- 
creases with decreasing energy. At the same time the 
ratio of velocities v±Tr/v\\o approaches one. The critical 
fixed point is, thus, characterized by an enhanced SU(2) 
symmetry. However, a peculiarity of the RG flow is that 
the velocity v«o itself vanishes in the low-energy limit re- 
sulting, e.g., in a diverging specific heat coefficient at the 
critical point. (In the opposite limit, v± 1l > v»q, run- 
away RG flow was found.) 

In the following, we address the question whether these 
critical properties are modified in the presence of a cou- 
pling to the spin degrees of freedom. 



B. Coupling to the spin sector 

Spin interactions in the one-dimensional Wigner crys- 
tal are described by the anti-ferromagnetic Heisenberg 
model with nearest-neighbor interactions. In the zigzag 
Wigner crystal, next-nearest neighbor interactions as well 
as ring exchange processes become important and lead 
to rich spin physics^ However, these additional interac- 
tions become important only once the lateral extent of 
the crystal is sufficiently large. Close to the transition, 
they are negligible, and the spin interactions are still de- 
scribed by the Heisenberg Hamiltonian, 



111) 



where the coupling constant J is exponentially small in 

Due to the exponential dependence of the spin inter- 
actions on the inter-particle distance, fluctuations of the 
electron positions immediately result in a modulation 
of J and, therefore, give rise to a magnetoelastic cou- 
pling between the spins and the charge modes. Specif- 
ically, the interaction energy Jj between electron j and 
j + 1 depends on the position of both electrons^ Jj = 
J(rj,r J+ i) ~ Jflrj+i — r 3 -|). In order to investigate the 
coupling between spin and charge modes perturbativcly, 
we expand the interaction energy Jj in small fluctua- 
tions around the equilibrium positions of the electrons us- 
ing 

Thus, 



ig \ r i+i- r j \ = gV(1 + u xj+1 - u xj ) 2 + (u yj+ i - u yj ) 2 . 



Jj ~ J(a) + aJ'(a) ( u 



xj+r 



if 



Note that the expansion of the coupling J in the longi- 
tudinal fluctuations u X j starts in linear order, yielding 



T~L 1 



-9\\ ^( u xj+i - u XJ )Sj ■ Sj+i, 



(12) 



with = — aJ'(a). By contrast, due to the symmetry of 
the one-dimensional Wigner crystal, the expansion in the 
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transverse fluctuations u y j begins only in second order, 
i.e., 

Hi = -g ± J2(u yj+1 - u yj ) 2 S, ■ S j+1 , (13) 



with g± = —aJ'(a)/2. 

The linear coupling to the longitudinal mode (fT2|) is 
familiar from the spin-Peierls problem^ In particular, 
the mode u\ with momentum q = w, u X j = u\{— 
couples to the staggered part of Sj ■ Sj+±, 

w« c «-2 ff |i5^ u |(-iy4-.^ +1 . (14) 



If the ul mode was sufficiently soft, this term would lead 
to a spin-Peierls transition. The crystal distorts such 
that the mode u\ assumes a non-vanishing expectation 
value giving rise to an alternation of weak and strong 
bonds, J ± 5 J. The system then gains magnetic energy 
by forming singlets on the strong bonds.— In our case, 
however, the magnetic energy is exponentially small such 
that it never can compete with the charge gap of the u\ 
mode, that is on the order of (n/m)y 1 nj 'as ln(n<i). We 
can, thus, conclude that the interaction of the spin 
degrees of freedom with the longitudinal modes does not 
influence the critical properties of the charge sector. 

We now turn to the coupling to the transverse mode 
(|13p . The most singular contribution is attributed to the 
critical mode. Substituting u y j = u^(— l) 3 , we obtain 



(15) 



It turns out, however, that this interaction is also irrele- 
vant as far as the critical properties of the Ising transition 
is concerned. This conclusion follows from a straight- 
forward power counting analysis of the Ising operator, 
(u^) 2 ~ 'J^VE', and the non-staggered spin-spin operator, 
Sj ■ S j+ i. 

Both magnetoelastic couplings, gu and g±, therefore 
do not modify the transition to the zigzag Wigner crys- 
tal described by the model ([8]). The transition happens 
only in the charge sector and the spin sector acts as a 
spectator. The spin and charge degrees thus remain es- 
sentially decoupled at the transition and the analysis of 
Ref. |3l| still applies. Wc now turn to the analysis of the 
transition in the weak-coupling limit nas 3> 1. 



IV. MULTIBAND QUANTUM WIRE 

At large densities, nas > 1, there exists a regime at 
sufficiently small oscillator length, xq < \J as /n, where 
the quantum wire contains well-defined single-particle 
subbands, see Fig. [2] In this regime, we can use the 
single-particle basis of the Hamiltonian (|T|), deriving from 



the kinetic energy and the confining potential, T + V^ on f, 
and expand the interaction Vint i n that basis. 

In second quantized form, the resulting Hamiltonian 
reads 



"-Els;-* 



nkcr nka 



(16) 



£ 

Wl,W2,Tl3,n4 

k,k' ,q;cr,(T / 



i(?) C n! 



k+qa n2h' —qa' n^k 1 a' n^ka ' 



where the electron operators c n ka destroy an electron 
with subband index n (i.e., the quantum number of the 
harmonic oscillator defined by V con f), momentum k in re- 
direction along the wire, and spin a =t, I- Two consecu- 
tive chemical potentials differ by the oscillator frequency, 

The electrons interact with an interaction amplitude U 
that depends on the transfered (longitudinal) momentum 
q and the subband indices rii . Its value is given by matrix 
elements of the screened Coulomb interaction in the basis 
of oscillator wavefunctions, 

/dq 
-^U{q,q y )T («»)> (17) 

where the Fourier transform of the interaction is given 
by 

U(q,q y ) = J dve-^U{v) (18) 

e y/q 2 + 1y 2 V ' 

and the matrix elements read 

^n 1 n 2 n 3 n 4 {%) = (19) 

dyi dy 2 eMy^^ ni { Vl W n2 {y 2 ) <j> n3 (y 2 )</>„ 4 ( Vl ) 



with the n th ID oscillator wavefunctions 4> n {y). Restrict- 
ing ourselves to the low-energy properties of the system, 
only interaction matrix elements U ni n 2 n 3 n i where the in- 
dices rii are pair-wise equal will appear. 

In the following, we focus on the situation where the 
first subband is filled, \x\ > 0, and the density of electrons 
in the second subband is very dilute, \fj, 2 \ <C fJ-i, i.e., the 
second subband is close to the quantum phase transition 
occurring at /i 2 = 0, see Fig. HJa). In this regime, we 
can neglect all higher subbands so that the band index 
is restricted to n = 1,2. The corresponding oscillator 
eigenfunctions are 



1/4 



1/4 



<fa(y) = V2mfi ( j ye 



(20a) 
(20b) 
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The reduced Hamiltonian for the two subband system is, 
thus, 



ft — Hi + H2 + rtl2) 



(21) 



where Hi represent the two (interacting) subbands, i = 
1,2, and H12 captures the inter-subband interactions. 
The condition x < \/ a B/n, that is fulfilled below the 
dotted line in Fig. [2J implies that the dimensionlcss 
interaction viU nin2n3ni {q) < 1 is small, where V\ = 
(2/ir)y/m/ (2/j.i) is the density of states of the filled first 
subband. 



A. Low-energy limit of the Hamiltonian 

In the following, we discuss the different parts of the 
two-subband Hamiltonian (|21[) separately before turning 
to the analysis of the full Hamiltonian in subsequent sec- 
tions. 



1. First subband 

The Hamiltonian Hi in Eq. (f2~Tj) represents the elec- 
trons in the first subband at finite [i\ > 0, interacting via 
the potential Unn(q). At low temperatures T -C p,\, we 
can approximate the first subband by a Luttinger liquid, 
linearizing the dispersion around the two Fermi points 
at ±kpi = ±y / 2m[ii. Introducing right- and left-moving 
fields, denoted R and L, respectively, 



c la (x) 



Jk F1 x 



—ikpix 



Li*(x), (22) 



the Hamiltonian Hi takes the standard form 



Hi = dx 



- iv P1 (rL&Ru, - L ld* L i*) (23) 
gicPmPiL — gisSiR.SiL ■ 



Here the particle and spin densities of right (r = R) and 
left (r = L) moving electrons are given by 



Ph 



= J2 r L r io, ^ir =\^r\ a a IJ ^r lal , (24) 



where Oi are Pauli matrices. The interaction constants 
describing scattering processes in the vicinity of the 
Fermi surface are given by 



9ic = C^ini (0) 



U llu {2k F1 ) 



g u = 2£/ mi (2fe F1 ) 



(25) 



Their magnitude is evaluated and discussed in 
Sec. IIV A 41 Here chiral interactions which only renor- 
malize the velocities of the charge and spin modes are ne- 
glected. The Hamiltonian Hi describes decoupled charge 



and spin modes, and has been intensively studied*^ The 
charge mode is gapless whereas the fate of the spin mode 
depends on the sign of gi s . If gi s were negative, the spin 
mode would acquire a gap. In the present case, however, 
gi s is positive, see Eq. (f25|) . and, therefore, for a single 
subband quantum wire, the spin sector remains gapless. 



2. Second subband 

The electrons in the second subband are close to the 
band bottom, [ii sa 0, and a clear distinction of left- and 
right-movers and subsequent linearization of the spec- 
trum is in general not possible. Furthermore, due to 
the low particle-density at the band bottom only inter- 
actions between electrons of different spin polarizations 
are of importance. Local interactions among electrons 
with the same spin are suppressed due to the Pauli prin- 
ciple. The residual interaction contains gradients of the 
electron fields; it is irrelevant in the RG sense at the crit- 
ical point fj,2 = and, therefore, will be neglected. The 
resulting Hamiltonian for the second subband reads 



H 2 = J dx[j2cl(x) (-JI-M2) c 2a (x) (26) 

+ v 41-0)44. i x ) c n (x)c n {x) , 



where 



V = 1/2222(0). 



(27) 



For small /Z2, the Hamiltonian Hi represents a strongly 
interacting system even for infinitcsimally small interac- 
tion V. This is best seen for negative chemical potential, 
/i2 < 0, where the ground state of Eq. (|26| is empty of 
particles. In this case, the retarded two-particle Green 
function, 

V 2 {x-x',t-t') = (28) 

- i9(t - t')([c n {x, t)c n {x, t), cl^x', t')c\ t (x\ t')]), 

can be straightforwardly evaluated at T = because it 
only requires the solution of a two-particle problem. The 
bare retarded two-particle function, T>^ , in the absence 
of interactions has the form 



+ i0 



2 V e + iO' 
(29) 



-2/i2 is the distance to the two-particle 



where e = w-^- 
mass shell. 

In the presence of the interaction V, the two particles 
repeatedly scatter off each other. The exact retarded 
two-particle Green function is then obtained by summing 
the interaction ladder in the particle-particle channel, see 
Fig. 02 resulting in 



_ T)(0)-1 



(k,u)-V. 



(30) 
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(a) 



T 



T 



(b) 

FIG. 3: (a) Two-particle Green function (|30[) of two elec- 
trons with opposite spin polarizations in the second subband. 
(b) The repeated scattering of the two electrons leads to the 
T- matrix, Eq. (|31[) . where the wiggly line represents the in- 
teraction V. 



Alternatively, this Green function can be expressed 
in terms of the two-particle T-matrix, T>2(k,ui) = 



V { ° } (k,uj)+V ( 2>(k,uj)T(k,uj)V^>(k,Lj), where 



,(o) 



2 0) ' 



T(k,u) 



V 



(31) 



Due to the inverse square-root singularity of T> 2 °\k,Lo), 
the T-matrix takes a universal form as the two-particle 
mass shell is approached, e — > 0, 



T(k,Lo) 



-2i 



SO 



(32) 



This corresponds to unitary scattering with phase shift 
S = 7r/2. In particular, note that the same limit obtains 
from Eq. pip for infinitely strong repulsion V — > oo. As 
the low-energy limit k, lo — > close to criticality fj,2 = 
coincides with the limit e — > 0, the second subband is 
populated by a strongly interacting electron gas (irre- 
spective of the value of V), where two electrons with 
opposite spin polarization cannot occupy the same state. 
As a consequence, the two-particle wavcf unction has not 
only nodes for electrons with the same spin as required 
by the Pauli principle, but also for electrons with op- 
posite spin polarizations. Such a strongly interacting 
gas has been dubbed an impenetrable electron gas£2-~— 
The impenetrable electron gas is the fixed point theory 
that describes the quantum phase transition in the sec- 
ond subband at /i2 = in the absence of inter-subband 
interactions. 

At positive /x 2 , the ground state of the system is a 
Luttinger liquid. For small < /12 < E p = mV 2 , the 
strong correlations of the emerging impenetrability of 
electrons are reflected in largely different energy scales 
governing the charge and spin sector. While the charac- 
teristic energy scale for charge excitations is the Fermi 
energy /12, the one for spin excitations, J, arises from 
corrections to the unitary limit (j3"2")l and is given by 3 - 3 . 
J ~ *C P2, where z/2 is the density of states 

in the second subband. As a consequence, the system 
remains in the Luttinger liquid state only at tempera- 
tures T < J. In the intermediate temperature range, 
J < T < fi2, spin excitations are incoherent, and the sys- 
tem can be described again as an impenetrable electron 




FIG. 4: Pair-tunneling process between the two subbands. 



gas. This phase has also been dubbed a "spin-incoherent 
Luttinger liquid" in the literature ! 33 ' 34 ' 63 

Thus, the Hamiltonian (|26|) close to the quantum phase 
transition at \ii = is governed by the physics of the 
impenetrable electron gas almost in the whole phase dia- 
gram except in the spin-coherent regime, T < ^2/(^2^), 
where the large spin entropy is released. 



3. Interaction between subbands 

Interactions between the subbands can be distin- 
guished into three types. There is a repulsive density- 
density interaction of electrons between the subbands 
with amplitude u c > 0, and a spin density-density in- 
teraction that is generically ferromagnetic, u s > 0. In 
addition, the amplitude ut describes pair-tunneling pro- 
cesses between the two subbands, sec Fig. @] The corre- 
sponding Hamiltonian rL\2 reads 



E \u c p 2 {RIB. + Plh) - U S S 2 ■ (SlR + SlL 

+u t J {c2 a cl s Li s Ri a + h.c.) (33) 



with 



-a. Here p 2 = c L c 2<r and S 2 = 
analogous to Eq. (f2~4")l . The values 



2 Eo-,er' c 2a^^^' C ^' ' 

for the couplings in terms of the interaction function U, 
Eq. (|17[) , are given by 

Uc = ^ r i22l(0)-it r i212(fcFl), U s = 2£/i212(&Fl), (34a) 

u t = U U 22{k F1 ), (34b) 

and are evaluated in the next section. As in Eq. (|26p . 
we neglected in (|33[) interaction processes involving ad- 
ditional spatial gradient terms. In particular, we disre- 
garded tunneling of electron pairs with the same spin 
polarization. 

The analysis of the full Hamiltonian (|21[) is compli- 
cated by the fact that the quantum critical point at 
fj,2 = for the interacting system has multiple dynamical 
scalesi 64 ' 65 Whereas the Luttinger liquid has a dynamical 
exponent z = 1, the spectrum of the second subband is 
characterized by z = 2. The multiple scales lead to the 
appearance of two different types of infra-red divergences 
in perturbation theory. The linear spectrum of the first 
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e(k) 



© ® 




FIG. 5: Energy spectrum as a function of the longitudinal 
momentum k of a two subband quantum wire with a lower 
subband 1 and a higher subband 2. We consider the quan- 
tum phase transition when the chemical potential reaches the 
bottom of the second subband, ^2 = 0. Depending on the po- 
sition of chemical potential 112 , we apply different approaches 
in sections ITVCl ITVDl ITVEl ITVFl and ITVGl as indicated in 
the figure. The energy scale E p is defined in Eq. (1391 1, 



subband, z = 1, yields logarithmic divergences while the 
quadratic spectrum of the second subband, z = 2, yields 
square-root singularities. In order to understand the in- 
terplay between these two types of divergences, we start 
in Sees . II V 51 II V D I by considering the case fj,2 < 0, where 
the two-particle Green function of the decoupled second 
(empty) subband is known exactly, see Eq. (|30p . Then, 
in Sec. IIV El we discuss the behavior of the system at 
the transition, fi2 ~ 0, and finally address the regime 
[ii > in Sec. IIV Fl and Sec. IIVGI The various regimes 
are indicated in Fig. [SJ 



4- Magnitude of the coupling constants 

Before turning to the analysis of the full effective low- 
energy Hamiltonian, let us determine the magnitude of 
the various coupling constants close to the phase tran- 
sition. The Hamiltonian (|2ip contains six coupling con- 
stants. For nas > 1, the transition happens at uxq ~ 1, 
i.e., kpl ~ .To. Both these length scales are smaller than 
the distance to the gate, kj,\ ~ i « d. Due to the long- 
range nature of interactions, the interaction constants 
corresponding to processes with zero momentum trans- 
fer are enhanced by a large logarithm. Furthermore, the 



dependence of C/„ iri2n3ri4 on the subband indices becomes 
negligible as d ^> xq- Thus, 



g ls ,u s ,u t 



gic V, u c 



v\naB 



^InA 



v\nas xq 



(35a) 
(35b) 



Note that, when the Wigner crystal regime is ap- 
proached, nag — > 1, the dimcnsionlcss couplings 
fig become of order one (apart from the logarithmic 
enhancement^), as expected. 



B. Perturbative analysis in the dilute limit /12 < 

If the chemical potential of the second subband is neg- 
ative, /i2 < 0, the ground state in the absence of inter- 
subband interactions consists of a Luttinger liquid in the 
first subband, Eq. (|23[) . and an empty second subband. 
In the following, we consider the perturbative renormal- 
ization of the intra- and inter-subband interactions. The 
perturbation theory in the inter-subband interactions is 
particularly simple because particle-hole polarizations of 
the second subband vanish exactly at T = 0. This allows 
to study the approach to the quantum phase transition 
for \i2 — > 0~. 

Perturbation theory encounters various singular cor- 
rections that diverge in the low-energy limit for [ii — > 0~ . 
Below, we list these singular correction to the various 
couplings given at some finite energy scale w. As men- 
tioned above, due to the multiple dynamics of the criti- 
cal point, different types of divergences, logarithmic and 
square-root singularities, can be distinguished. The type 
of divergence depends on the intermediate state. A 
square-root divergence requires that in the intermediate 
state both particles are in the second subband. If at least 
one particle in the intermediate state is in the first sub- 
band, a logarithmic divergence is obtained. Here one can 
further distinguish between logarithmic divergences that 
are cut off by the chemical potential [12 , if one of the par- 
ticles in the intermediate state is in the second subband, 
and those that are not, when both particles are in the 
first subband. 

The perturbative corrections to the intra-subband in- 
teractions in the occupied first subband are given by 



f 



he 



5gi 



9ls 



1 ^0 

■ in r 

2itvfi M 



2w?r4 0) (0,u;) 



(36a) 
(36b) 



where Eq is a UV cutoff on the order of the Fermi en- 
ergy of the occupied first subband, Eq ~ Ep\. Apart 
from the standard logarithmic corrections^ there are 
additional singular corrections due to the pair-tunneling 
Ut between the subbands. The corresponding diagram 
is shown in Fig. [BJa). These singular corrections yield 
square-root divergences as they involve the bare two- 
particle Green function V^(0, uj), see Eq. (p?9"|) . that be- 
haves as 1/^/^7+27*2 • 

The interactions of particles within the empty subband 
obtains the singular correction 



SV 



u 2 t , E 
In—, 



(37) 



The first term on the right hand side corresponds to the 
first term in the expansion of the two-particle T-matrix 
(f3"TT) in powers of V, see Fig.[B][b). The second term due 
to pair-tunneling to the first subband lowers the effec- 
tive interaction. It is represented by a diagram similar 
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(a) 



(b) 



FIG. 6: Vertex corrections that involve the two-particle Green 
function of the second subband and diverge as 1/s/e with 
vanishing energy e: correction to (a) the intra-subband in- 
teractions g\ x , x = c, s; (b) the intra-subband interaction V, 
and (c) the pair-tunneling u t . The dashed/solid lines repre- 
sents the propagator in the first /second subband, the dot is 
the pair-tunneling amplitude u t , and the wiggly line is the 
intra-subband interaction V. 



to Fig. [UJa), but with the dashed and solid lines inter- 
changed. 

Finally, the corrections to the inter-subband interac- 
tions read 



Su c 
8u.« 



■ In ■ 



E 



2-KVFl \w + fJ,2\' 



In 



2ttv F i \u + fi 2 \' 



(38a) 
(38b) 



Su t =Vu t V^>(0,Lo) 



,(0) 

2utu, 

TTVFl \u + H2\ 



2txvf\ 



hi- 



In. 



En 



The two-particle Green function of the second subband 
enters also here, namely in the correction to the pair- 
tunneling vertex, see Fig.[5][c). 

The singular vertex corrections listed above are either 
regularized by a finite energy w or a finite chemical po- 
tential ^2- We will first focus on the singular corrections 
that survive in the limit of large negative \i2- 



C. Perturbative regime at large negative ii 2 - dilute 
weakly interacting Fermi gas in the second subband 

For large negative chemical potential of the second sub- 
band, the vertex corrections involving the two-particle 
Green function of Fig. [5] can be treated perturbativcly. 
The relevant energy scale is found by estimating the 
value of the square-root corrections in the limit uj — > 
and requiring Sx/x < 1, where x = g lc , g ls , V, Ut- 
Using 2? 2 °^ oc ^/m/yu 2 , one obtains \^\ 3> E p = 
max{mMj 1 9i c , mu\ /g 2 s , mV 2 }. Thus, for our model pa- 
rameters, Eqs. (|33|. the energy scale E p is determined by 
the intra-subband scattering V, 



E v = mV' 



(39) 



Similarly, for |^2| 3> E p , the logarithmic corrections in- 
volving the chemical potential fj,2 in their argument will 
reduce to small perturbative corrections only. The re- 
maining logarithmic corrections to gi s , V, u t that are sin- 
gular in the limit w — > derive from exciting the Lut- 
tinger liquid and can be summed up with the help of the 



conventional RG approach. After integrating out modes 
within an energy shell (E /b, E ) with the scaling param- 
eter b > 1, we obtain the RG equations 



dg 



is 



ah 



dhib 
dV 

dlnb 
du t 

<91n6 



2'KVpi 



TTVFl 

(gic + 3gi s /4)u t 

2'KVpi 



(40a) 
(40b) 
(40c) 



The solution of these equations is straightforward. The 
first equation is the standard RG equation for the 
spin mode of a Luttinger liquidJ^ The coupling g ls is 
marginally irrelevant so that it vanishes logarithmically, 
gis{b)/ (ttvfi) (x l/ln&, for vanishing energies, b — >• oo. 
The last two equations describe the renormalizations of 
the interaction V in the second subband and the pair- 
tunneling ut due to the interaction with the Luttinger 
liquid. 

Although the latter two interactions, V and Ut, involve 
high-energy excitations to the second subband that is far 
away in energy (\fi2\ 3> E p ), it is nevertheless impor- 
tant to consider them in order to check for consistency. 
In particular, if |V| flows to large values, the condition 
(38c) \n2 1 3> E p = mV 2 may be violated during the RG flow. 

The pair-tunneling Ut leads to a reduction of the inter- 
action V and could potentially drive it to large negative 
values. A change of sign of V from repulsive to attrac- 
tive is accompanied by the appearance of a two-particle 
bound state at energies below the bottom of the second 
subband. In that case, the transition at /i 2 — could be 
preempted by the condensation of bound pairs. 

However, it turns out that the RG flow of the interac- 
tion V is rather short. The pair-tunneling u t which drives 
the flow of V suffers an orthogonality catastrophe, i.e., it 
is suppressed by the interactions in the Luttinger liquid. 
Thus, Ut is irrelevant and flows to zero with a scaling di- 
mension that approaches g\ c j '(2ttvfi) in the low-energy 
limit, b — > oo. It is straightforward to solve Eqs. (|4"0"|) for 
the effective coupling V e g at the lowest energies. Namely, 



gic 



> 9u 

"9u 



(41) 



where the function <j> has the limits 0(0) = 1 and <fi(x) 
4/x for x — > oo. Using the initial values of the interaction 
constants, Eq. (|35|) . one obtains SV/V ~ l/(ln<i/xo) 2 , 
i.e., for d ^s> x the correction is logarithmically small. 
Thus, the effective interaction V e g remains repulsive and 
in the perturbative regime, and, consequently, the condi- 
tion \fi 2 \ m Kff s ^i n holds. 

To summarize, in the limit of large negative chemical 
potential, |/l*2 | 3> E p , the ground state of the system is 
not affected by the inter-subband interactions which re- 
main perturbative. The first subband is a Luttinger liq- 
uid while the second subband corresponds to a weakly- 
interacting dilute Fermi gas whose population is expo- 
nentially small in ^/T for | /X2 1 3> T. 
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D. Quantum phase transition at /U2 = 0: Lifshitz 
transition of impenetrable polarons 



If the negative chemical potential ^2 < increases and 
passes the threshold [/X2 1 ~ E p , see Eq. (|39[) . the square- 
root singularities of the two-particle Green function T> 2 °^ 
develop in the perturbative corrections to the couplings 
gic gis, V, and u t , see Eqs. (|36l37l38c|) and Fig. [6j In the 
following, we consider the limit of small negative chemical 
potential, | 1 ^ E pi and, in particular, the approach to 
the quantum phase transition fi 2 — . 



1. Formation of an impenetrable electron gas in the second 
subband 



(a)-' 




(b) - 

FIG. 7: (a) Renormalization of the pair-tunneling vertex ut- 
Two electrons of the first subband (dashed lines) are trans- 
fered to the second subband (solid lines). In addition to the 
bare pair-tunneling (black dot), the vertex receives renormal- 
izations from the repeated scattering of electrons in the second 
subband that is captured by the T-matrix, see Eq. (|31(l and 
Fig. [3] (b) Renormalization of the intra-subband interactions 
gi x with x = c, s within the first subband. The black square 
represents the bare interaction; the renormalization is due to 
pair-tunneling into the second subband with the renormalized 
pair-tunneling vertex defined in (a). 



In the limit \fi 2 \ <C E p , the inverse square-root singu- 
larities that appear in the perturbation theory become 
of order one for energies lo ~ E p , while the logarithmic 
corrections are still small. So it is permissible to consider 
first the stronger inverse square-root singularities deriv- 
ing from the two-particle Green function of the second 
subband. In order to identify the class of most diver- 
gent diagrams, we have to consider processes that con- 
tain as many two-particle propagators T>2 as possible. 
These processes are associated with the repeated scatter- 
ing of particles in the second subband. The summation 
of this class of diagrams corresponds to the formation of 
an impenetrable electrons gas in the second subband as 
discussed in Sec. IIV A 21 As a consequence the bare two- 
particle propagator T>^ appearing in Eqs. (|36l37l38c[) 
has to be replaced by the full two-particle propagator 

Consider first the renormalization of the intra-subband 
interaction V in the second subband. The summa- 
tion of repeated two-particle scattering processes pro- 
motes the interaction V to the T-matrix, V — > V e g = 
V + V 2 V 2 (k,uj) = T(k ,u>), see Eq. (3T]) and Fig. [3] As 
discussed in Sec. IIV A~2l the effective interaction thus be- 
comes explicitly energy dependent and, as the distance 
e = oj — k 2 / (Am) + 2fi 2 to the two-particle mass shell fur- 
ther decreases, |e| < E p , it assumes the universal form 

T(k,uj) w — 1/T>2 which is characteristic for an 

impenetrable electron gas. 

The formation of an impenetrable electron gas in 
the second subband has the consequence that the pair- 
tunneling becomes inefficient because the probability to 
find two electrons at the same position in space is strongly 
suppressed. The singular vertex corrections for the pair- 
tunneling are shown in Fig. [TJa). In the limit e — > 0, 
summation of these diagrams gives 

u t (e) = u t + Vu t V 2 (k,cj)^ ^ . (42) 



Thus, the effective pair-tunneling vanishes as 



u t (e) ~ u t \ — fro for e -> 0, (43) 
V mr 

i.e., pair-tunneling is suppressed by a factor e / E p . 

Furthermore, the renormalizations of the Luttinger liq- 
uid interactions due to pair-tunneling, Eqs. (|36|) . are now 
regularized by the strong interactions in the second sub- 
band, see Fig.Efb). Replacing vf ] by V 2 in Eqs. (36]) 
and taking the limit e — > 0, one obtains 

gf x =9i x -C x y, (44) 

where x — c, s and the constants C c — 1/2, C s = 2. 

Thus, the pair-tunneling ut leads to a negative cor- 
rection to the Luttinger liquid interactions. If u t is suf- 
ficiently large, it may lead to a sign change of the in- 
teraction constants and, thus, induce attractive (anti- 
ferromagnetic) interactions in the Luttinger liquid. In 
particular, the smaller effective coupling gf^ , that con- 
trols the spin sector, turns negative if 

vgu 

If the effective coupling were driven negative, gjf < 0, 
the interactions between the residual low-energy degrees 
of freedom would drive the interaction to strong coupling 
and generate a spin gap in the first subband^ see also 
discussion in the next section. However, for the parame- 
ters of Eqs. (|35p we have u 2 /(Vgi s ) ~ l/]n(d/xo) so that 
the criterion (|45|) is never fulfilled for d xq. Thus, the 
Luttinger liquid in the first subband remains gapless. 

So far we considered the effect of repeated scattering 
of particles in the second subband yielding the most sin- 
gular corrections. We found that, as a consequence, an 
impenetrable electron gas forms in the second subband at 
energies u> ~ E p , and the pair-tunneling between the sub- 
bands becomes ineffective due to the formation of nodes 
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in the two-particle wavefunctions. However, the remain- 
ing interaction still leads to interesting physics that de- 
velops for even smaller energies, w < £ p , as we show 
below. 



2. Polaron formation at lowest energies 

Due to the formation of an impenetrable electron gas in 
the second subband, the pair-tunneling is suppressed by 
a factor yJuT/ 'E p , see Eq. (j43]). and consequently becomes 
ineffective for energies ui < E p . At the lowest energies, 
lu <C E p , the pair-tunneling may therefore be neglected. 
However, logarithmic divergences in Eqs. (|36b|38b| sur- 
vive and lead to a residual flow of the remaining vertices. 
A perturbative treatment of the vertex corrections at en- 
ergies oj < E p yields the one-loop RG equations 



(46a) 
(46b) 



dgis 


9u 


d\nb 


2itvfi 


du s 


< 


d\nb 


2-kvfi 



Whereas the first equation derives from processes in the 
Luttinger liquid of the first subband only, the flow of u s 
is caused by the virtual excitation of a particle to the 
empty second subband. Contrary to the pair-tunneling 
process, there is only a single particle involved which - as 
the second subband is empty - cannot repeatedly scatter 
from other particles. As a result, u s is not suppressed 
by the formation of an impenetrable electron gas in the 
second subband even for uj < E p . 

To solve the RG equations, it is sufficient to consider 
the logarithmic renormalizations arising at the lowest en- 
ergies, i.e., the effective cut-off can be chosen as E < E p . 
Consequently, the initial values for the RG flow are given 
by the values of the coupling constants at the energy scale 
Ep, namely the effective coupling gf^ of Eq. (|4"4"]l and the 
bare inter-subband coupling u s . 

Eq. (|46a[) is the standard RG equation for the Lut- 
tinger liquid. Depending on the sign of gf^ , the flow is 
cither to weak or strong coupling. As discussed in the 
context of Eq. (|4"4"|). for d/xo ^> 1 the coupling gff is 
positive, corresponding to ferromagnetic spin-spin inter- 
actions, and thus flows to weak coupling. 

Eq. (|46b[) for the inter-subband spin-spin interaction 
has the same form as Eq. ()46a|) . Thus, the (positive) cou- 
pling constant w s (w) also decreases logarithmically with 
decreasing w, 



u s (uj) 



1 



■hi 



(47) 
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i.e., the ferromagnetic inter-subband spin-spin interac- 
tion is marginally irrelevant. Contrary to the intra- 
subband spin-spin interaction, here the flow is stopped 
for a finite (negative) chemical potential |^t 2 | <C E p . 

Even though the pair-tunneling is suppressed by the 
strong interactions in the second subband, the low-energy 



FIG. 8: Self-energy for electrons in the second subband. The 
dot represents either the charge or spin inter-subband inter- 
action vertex, u c or u s , respectively. 



limit of the system does not correspond to two decou- 
pled subbands, namely a Luttinger liquid and an im- 
penetrable electron gas. The inter-subband charge and 
spin density interactions, u c and u s , respectively, lead 
to a polaron effect^ The particles of the second sub- 
band occupy states close to the band bottom, see Fig. [SJ 
and as a consequence, they propagate very slowly: their 
velocity vanishes close to the band bottom. By com- 
parison, the density fluctuations in the filled first sub- 
band are rather fast, and the density can adiabati- 
cally adjust to screen the particles in the second sub- 
band. This adjustment is reflected in logarithmically 
singular corrections to the residue Z of the fermionic 
single-particle Green function in the second subband, 
Q2 (u>, k) = cj— fc 2 /(2m)+fi 2 - S(w, k), that only appears 
in two-loop order. Evaluating the self-energy diagram 
shown in Fig. [51 one obtains 



S(w,fc) = \u) - - h^ 2 



2m 



2ui 



^v 2 



x In ■ 



(2^fi) 2 
E(j 



(48) 



OJ 



Note that Fig. [8J represents the only important two-loop 
self-energy diagram for /i2 < and T = 0. In particular, 
contributions arising from intra-subband scattering V are 
exponentially suppressed in H2/T, see Ref. [6l]- The log- 
arithmically singular self-energy £ results in an RG flow 
for the residue Z giving rise to an anomalous dimension 
for the electrons in the empty second subband. The flow 
can be integrated after also taking into account the two- 
loop vertex corrections^ Note that not only charge but 
also spin excitations (absent in the spin polarized case 
considered in Ref. |28| ) contribute to screening and conse- 
quently to the logarithmic singularity in Eq. (|48l) . 



3. Universality class of the quantum phase transition 

After all, we can identify the quantum phase transition 
occurring at [ii = as a Lifshitz transition corresponding 
to the filling of an empty subband as a function of the 
chemical potential ^2- The electrons that fill this empty 
subband are, however, strongly interacting, characterized 
by unitary scattering. Furthermore, each of these elec- 
trons is screened by the density and spin excitations of 
the Luttinger liquid in the filled first subband. As a re- 
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suit, the quantum phase transition is a Lifshitz transition 
of impenetrable polarons. 

In the following, we discuss a rcgularization scheme for 
the square-root singularities that allows for an unbiased 
renormalization group analysis and confirms the physi- 
cal picture for the quantum phase transition developed 
above. 



Regularization of the z — 2 singularities at the 
quantum critical point \i% = 



E. 



The two subband Hamiltonian, Eq. (|21j) . has previ- 
ously been considered by Balents and Fisher in Ref. [3fl 
in the context of the two-chain Hubbard model. In order 
to deal with the logarithmic and square-root singulari- 
ties encountered in perturbation theory, see section llVBl 
they considered a generalized dispersion for the electrons 
in the second subband, e(fc) = \k\ 1+e v 1 ~ e /(2m) e , where 
v is an artificial parameter with the dimension of veloc- 
ity. The physical quadratic dispersion of the Hamiltonian 
(|2T))) is recovered for e = 1 whereas e = corresponds to 
a linear spectrum. Treating e as a small parameter the 
square-root singularities are regularized and RG equa- 
tions for the Hamiltonian close to criticality, [ii w 0, can 
be derived. In lowest order in e, one obtains^ 



(49a) 
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(49c) 
(49d) 
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Here the RG equations are derived under the assump- 
tion that all dimensionlcss coupling constants are much 
smaller than e. 

The parameter e appears in the equations for the intra- 
subband interaction in the second subband, V, and the 
pair tunneling, u t , i.e., the coupling constants that ac- 
quire corrections due to the interaction V. Both V and 
Ut initially have a positive scaling dimension, e and e/2, 
respectively, i.e., they are both relevant. However, the 
scaling dimension of V is larger and, thus, the intra- 
subband scattering V is the most relevant perturbation. 
Initially, V therefore grows exponentially as a function 
of In b. As a consequence, the intra-subband scattering 
quickly reaches an attractive fixed point V/(2irv) — > e. 
This fixed point can be identified with the impenetrable 
electron gas within the second subband. 



The pair-tunneling initially increases as well but with 
the smaller scaling dimension e/2. As V approaches 
the strong-coupling fixed-point, the scaling dimension of 
the pair tunneling changes from relevant to irrelevant as 
e/2 — V/(2nv) — > —e/2, corresponding to a suppression 
of pair-tunneling due to the formation of nodes in the 
two-particle wave function, see also Eq. (|4"3")l . As a con- 
sequence, the pair-tunneling quickly flows to zero and the 
remaining RG flow reproduces the one given in Eqs. (|46|) . 

The approximations of section TlVD II with the expres- 
sions (|4"4"]) for the effective couplings are exactly recovered 
if the initial exponential flow of V and Ut is reduced only 
to the first two terms on the right hand side of Eqs. (|49c|) 
and (|49fj) . Plugging the obtained result for ut(b) into 
Eqs. (|49ap and (|49b[) . one recovers Eq. (|4"4"|) with the same 
coefficients C x . 

It is interesting that the reduction of the effective inter- 
actions (|44p quantified by C x can be mainly attributed 
to an approximate unstable fixed-point of the RG flow 
(|4*9"|) . The interactions g\ c and g\ s are substantially 
renormalizcd when the pair-tunneling reaches its max- 
imal value. This occurs at scales b ~ b m when the flow 
of ut is approximately stationary, i.e., when the intra- 
subband scattering reaches half of the fixed point value, 
V(b m )/{2irv) = e/2. Close to this stationary point, the 
flow of u t can be approximated by 



u t {b) w uta 



■ ln"(b/6 m ) 



(50) 



The maximal value of the pair-tunneling, u tmax , is ob- 
tained from the approximate RG invariant of the ini- 
tial flow. Considering the coupled equations (|49c|) and 
(|49fj) . and keeping only the dominant first two terms in 
Eq. (|49f[) . one finds the approximate RG invariant I = 
(eV-V 2 /(2TTv))TrvFi/u 2 + 2\nu t . With V(b m )/{2-Kv) = 
e/2, the maximal pair tunneling then obtains u tmax « 
Uf\/ivve/(2V), i.e., the pair-tunneling close to b m is en- 
hanced by the large factor y / nve/(2V). 

Consequently, the pair-tunneling correction will dom- 
inate the renormalization of g\ s close to b m so that we 
can approximate its flow by dg\ s j '9 In 6 m —u 2 /{ttv) with 
u t (b) given by Eq. (|50[) . Integrating this equation one 
obtains an effective intra-subband interaction 



eff 

9u 



9u - VTTy 



(51) 



The stationary point of the RG flow thus accounts al- 
ready for a reduction with a prcfactor y/ir as compared 
to C s = 2 obtained in section llV D 11 Needless to say that 
we arrive at the same conclusion concerning the stability 
of the Luttinger liquid as in the previous section, see the 
criterion of Eq. (|45|) . 

This is to be contrasted with the results of Ref. |H, 
where a Hubbard ladder was considered. In that case, all 
interactions are of the same order, namely on the order 
of the Hubbard interaction U . Thus, for Hubbard initial 
conditions the criterion (|4"5"j) is fulfilled, and a spin gap 
can be generically expected in agreement with Ref. l35l . 
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For the quantum wire with screened Coulomb interac- 
tion, however, the system remains gapless. 



F. Dense limit of large positive chemical potential, 

H2 > Ep 

The square-root singularities of the perturbation the- 
ory at fi2 = can also be regularized by considering a 
finite positive chemical potential larger than the strong- 
coupling scale, fi2 > E p , defined in Eq. (f3U)) . For a 
positive chemical potential, the ground state contains 
a finite density of particles in the second subband. At 
lowest temperatures, T U2, we can then linearize 
the quadratic spectrum of the second subband and map 
its Hamiltonian 7^2, Eq. (|26|) . to a Tomonaga-Luttinger 
model. 

Analogous to the treatment of the first subband 
in Sec. IIV A 21 we introduce right- and left-movers, 
cf. Eq. (|2"2"]l . to find the standard Luttinger Hamiltonian 
% 2 m ; cf- Eqs. (|23l24p . with interaction constants g2 C , <72s- 
In the dense limit, Epi 3> U2 3> E p , the values of the 
coupling constants are 



32c 



V 
2"' 



92. s 



2V. 



(52) 



In terms of the new fields, the inter-subband interac- 
tions (132|) can be written as H 1 ^ = H^ ca + with 



%i2 CS = u c j dx {p\RP2L + P-zrPil) (53) 

-U a J dx (jSlR ■ S 2 L + S 2 R ■ Sil 

and 



U% = I dx 



Utc [T R T L + h.c.) - u ts [T R -T L + h.c 



(54) 



where T r = J2a r i<j r 2<? and T r = Da r L^,^' r 2a'/ 2 - For 
fi2 "C Ep\, the charge and spin pair-tunneling interac- 
tions are related to the couplings defined in Eq. ([3"3"| by 
u tc = Ut/2 and u ts = 2u t . 

The resulting Hamiltonian Hi +n 2 in + 'Hi2, with %i of 
Eq. (|23[) . was previously analyzed by Varma and Zawad- 
owski in Ref. |36| and subsequently by other authors . 35 ' 37 
Integrating out modes in an energy shell (E^/b, Eq) with 
b > 1 , where the cut-off here is on the order of the small 
chemical potential, Eq ~ U2 <C Epi, one obtains the 
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Note that this pcrturbativc approach necessarily breaks 
down as the phase transition U2 — > is approached be- 
cause the vanishing Fermi velocity vf2 — \/2nifi2 causes 
terms on the right hand side of the RG equations to di- 
verge. The perturbative treatment is controlled as long 
as all the terms on the right hand side of Eqs. ([55]) arc 
smaller than the coupling constants they renormalizc. 
For the parameters of our model, Eqs. (|35|) . the largest 
term is ~ g2s/ v F2 in Eq. (|55d|) . Thus, we obtain the con- 
dition g2s/vF2 *C 1 or, equivalently, [12 3> E p with the 
strong-coupling energy scale E p = mV 2 of Eq. ((55)) . 

Close to the transition U2 > E p , the initial RG flow 
is determined by the terms that are enhanced by factors 
of the inverse Fermi velocity, Vp 2 £& Keeping only such 
terms, the flow can be approximated as 
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(56b) 
(56c) 
(56d) 
(56e) 



This initial flow shares similarities with the one obtained 
at large negative chemical potential I/Z2I 3> E p in section 
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IIV C[ but with the role of the two subbands reversed: 
here only processes contribute where both particles in the 
intermediate state are in the second subband. The equa- 
tion for the intra-subband spin interaction g2 S decouples. 
As g2 S > 0, the flow is to weak coupling. The pair- 
tunneling reduces the intra-subband interactions of the 
first subband, g\ c and g\ Sl making them less repulsive. 
This rcnormalization is, however, limited because the 
pair-tunneling again suffers an orthogonality catastrophe 
- now due to the interactions in the second subband. 
Eqs. (|56dl56ep show that the pair-tunneling is irrelevant 
with scaling dimension 32c/(27TOf 2) and, thus, quickly 
decreases. Introducing linear combinations of coupling 
constants, Eqs. (|56d|) and (|56e[) can be recast in the form 
du ta /d\nb = -g2aU ta /(2nvF2), where x a = x c + ax s 
and a = —1/4,3/4. Using the identifications with the 
coupling constants of the Hamiltonian (J2TJ) , one finds 
the initial conditions Ut(3/4) = 2it t = 2£/ 112 2(fcFi) and 
Mt(_i/4) = 0, so that the latter does not flow. Solving 
consecutively the RG equations for </2s, Mt(3/4), and then 
gi s , the finite renormalizations of the intra-subband scat- 
tering <7i s during this process can be written as 

off 1 {utc+juuf Mf g 2s ^ AMA\ u t 
9u = 9n - r> ~ n — ) = 9u ~ 4<M4) — 

(57) 

with the same function <fi as in Eq. (|4ip. Here (f)(4) « 
0.34. 

As above, the value of gff determines the fate of the 
subsequent flow governed by terms that are not enhanced 
by the inverse velocity: While for positive gf^ the flow 
is towards weak coupling, run-away flow obtains for neg- 
ative gf^, signaling the opening of gaps. The coupling 
changes sign if u^/(Vgu) > 0(1) so that we reproduce 
again the stability criterion (|4"5j) . 

In principle, the corrections to the RG flow due to 
terms that are not enhanced by Vp2 could induce other 
instabilities. The pair-tunneling, that is irrelevant dur- 
ing the initial flow (|56[) , renormalizes the charge coupling 
32c and could drive it attractive. This could, in turn, 
make the pair-tunneling relevant, resulting in an insta- 
bility. However, in our case the criterion for this second 
scenario of run-away flow is even weaker than Eq. (|45|) . 
Thus, we conclude that in the regime [12 > E p the ground 
state of two Luttinger liquids in the two subbands is sta- 
ble. 



G. Small positive chemical potential: < fi2 < E p 

So far, we discussed all regimes in Fig. [5] close to the 
phase transition except the one for finite positive but 
small chemical potential, < fj,2 < E p . The physics in 
this regime is already strongly influenced by the strong- 
coupling fixed point of the impenetrable electron gas that 
controls the quantum phase transition. Whereas in the 
regime of small and negative ^2 of section IIVDI a theo- 



retical description in terms of only two interacting parti- 
cles within the second subband was sufficient, now a fi- 
nite density of strongly- interacting electrons populate the 
second subband which eludes a perturbative description. 
As the charge sector of the impenetrable electron gas 
corresponds to a spinless Fermi gas^ we can only spec- 
ulate that the physics within this strongly-interacting 
regime shares certain similarities with a system of spin- 
less fermions. 

In a two-subband system of spinless fermions, one 
generically finds an instability for small but finite H2 
because pair-tunneling becomes relevant and leads to 
the opening of a gap 3^ The main reason is that, for 
the spinless case, the intra-subband interaction of the 
weakly populated second subband is of order 0(v 2 F2 ) 
due to the Pauli principle, whereas pair-tunneling is only 
suppressed by a factor vf2, namely Ut = t/ii22(fcFi — 
&F2) — ^1122(^1 + kp2). As a result, the orthogonality 
catastrophe arising from adding particles to or removing 
particles from the second subband is substantially weak- 
ened: the RG flow does no longer contain terms of order 
Vp2 as in Eqs. (|55[) while the remaining terms make the 
pair-tunneling marginally relevant. The flow of the pair- 
tunneling to strong coupling signals the appearance of a 
gap due to a locking of the relative charge mode. Im- 
portantly, for spinless fermions there does not exist any 
energy scale E p that limits the range of validity of the 
perturbative RG in contrast to Eqs. (|5"5|) . 

It is not unlikely that a similar kind of relevant pair- 
tunneling mechanism might be at work in the spinfull 
case as well in the regime < [12 < E p where the physics 
is already influenced by the impenetrable electron gas 
fixed-point, maybe giving rise to a gap in the charge sec- 
tor. 



V. SUMMARY & DISCUSSION 

We analyzed the quantum phase transition in a 
quantum wire from a one-dimensional to a quasi-one- 
dimensional state assuming that electrons interact via a 
long-range Coulomb potential screened by a nearby gate. 
The model of the quantum wire was defined in section HI] 
and a mean-field phase diagram was presented in Fig. [5] 

In the limit of strong interactions, uclb *C 1, the 
ground state is determined by the charge degrees of free- 
dom. The charges minimize the Coulomb energy by form- 
ing a Wigner crystal. Exchange processes between elec- 
trons give rise to an exponentially small exchange cou- 
pling between the spins so that the spin sector is well 
approximated by the nearest-neighbor Hciscnbcrg model. 
The phonons of the Wigner crystal and the spin excita- 
tions of the Heisenberg model thus account for a gapless 
charge and spin mode in the one-dimensional limit giv- 
ing rise to a C1S1 phase in the notation of Ref. HH. At 
the transition to a quasi-one-dimcnsional state, the ID 
Wigner crystal deforms and splits into two rows, see sec- 
tion |TlT] and Fig. [T^b). The transition is characterized 
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by a (non-local) Ising order parameter. The exponen- 
tially small spin exchange between nearest neighbors de- 
pends sensitively on the distance between charges which 
results in a magnetoelastic coupling between the spin and 
phonon excitations of the Wigner crystal. We found that 
due to the Ising symmetry the spin exchange can only 
couple inefficiently to the square of the order parameter 
so that the universality class of the transition remains 
unaffected by the spin sector. As a consequence, the 
analysis of the spin-polarized case also applies herej^I 
and the transition is an Ising transition with logarith- 
mic corrections resulting from the coupling of the critical 
Ising modes to the longitudinal phonons. Although this 
latter coupling was found to be marginally irrelevant, its 
RG flow leads to an enhanced SU(2) symmetry accompa- 
nied by a suppression of the velocity of the excitations ^ 
Right after the transition to a deformed zigzag Wigner 
crystal, the resulting state possesses a single gapless lon- 
gitudinal phonon and a gapless spin mode, corresponding 
to a CIS 1 phase. If the electron density increases further, 
the fate of the spin sector depends on the strength of ring 
exchange processes^ 

In the limit of weak interactions, nas 3> 1, the transi- 
tion can be analyzed in terms of a two-subband model 
of interacting electrons, see section IIVI In the one- 
dimensional limit, only the lowest subband is filled, form- 
ing a Luttingcr liquid with a gapless charge and a gap- 
less spin mode, i.e., it is a C1S1 phase. As the density 
increases, a second subband starts being populated, see 
Fig.[TJa). This quantum phase transition is characterized 
by an interplay of charge and spin degrees of freedom and, 
more importantly, multiple dynamical scalesj 31 ' 64 ' 65 The 
Luttingcr liquid within the first filled subband is Lorentz- 
invariant and thus has a dynamical exponent z — 1. The 
electrons in the empty second subband, on the other 
hand, are characterized by Galilean invariance and the 
dynamical exponent is z = 2. As explicitly shown in 
section IIVB1 the presence of two dynamical exponents 
is reflected in divergences in perturbation theory of dif- 
ferent strengths. Whereas perturbative corrections in- 
volving excitations of the Luttinger liquid are accompa- 
nied by ubiquitous logarithms, the corrections involving 
electron-electron polarizations in the second subband are 
more singular and diverge as the inverse square-root of 
energy. Indeed, the electrons at the bottom of the sec- 
ond subband are strongly interacting. This is best un- 
derstood by considering the dimensionless coupling U2 V, 
which diverges at the band bottom due to the diverging 
one-dimensional density of states v<i. As a result, the 
scattering between two electrons of different spin polar- 
izations is unitary, and the electrons become impenetra- 
ble. We found that this emergent impenetrability in fact 
stabilizes the Luttinger liquid in the filled subband for the 
following reason. The system is able to gain energy by 
pair-tunneling processes of electrons with different spin 
polarizations from the filled first into the second empty 
subband. This effectively generates an attractive inter- 
action in the Luttingcr liquid that increases the tendency 



towards the formation of a spin gap. However, due to the 
emergent unitary scattering the probability to find two 
electrons at the same position in space decreases dramat- 
ically and pair-tunneling between the subbands becomes 
inefficient. We showed that for electrons interacting via a 
screened Coulomb potential the gapless spin excitations 
within the first subband survive at the transition, be- 
cause the bare pair-tunneling between the subbands is 
weaker than the intra-subband interaction, see the crite- 
rion in Eq. (|45|l . This is in contrast to the two-chain Hub- 
bard model where both couplings are of the same order 
and pair-tunneling promotes the formation of bound elec- 
tron pairs^ We showed in section lTV El that the stability 
criterion of Eq. P5|) , that embodies the competition be- 
tween intra-subband interaction and pair-tunneling, can 
be attributed to a stationary, i.e., a repulsive approxi- 
mate fixed-point of the renormalization group flow. In 
addition, we found similarly as in the spinless case that 
the impenetrable electrons are screened by plasmon and 
spinon excitations in the first subband. The quantum 
phase transition at weak interactions is, thus, identified 
as a Lifshitz transition in terms of impenetrable polarons. 

We were not able to identify the phase and, in partic- 
ular, the number of gapless modes for a small but finite 
filling of the second subband, because close to the Lif- 
shitz transition the physics is still strongly influenced by 
the strong coupling fixed-point of the impenetrable elec- 
tron gas that is beyond a perturbative treatment. Wc 
conjecture that - similarly to the spinless case^ - pair- 
tunneling processes might play here an important role as 
well, maybe reducing the number of gapless modes. We 
conjecture that pair-tunneling processes - similar as in 
the spinless case2£ - might play here an important role 
as well, maybe reducing the number of gapless modes. 
If the density is further increased and the filling of the 
second subband is sufficiently large, the interactions can 
again be treated perturbatively and the filled two sub- 
bands form two decoupled Luttinger liquids (though not 
in the original basis) with two gapless charge and spin 
modes, i.e., a C2S2 phase. 

Another open question concerns the connection be- 
tween the strong and the weak coupling limit. There 
should exist a multicritical point at intermediate den- 
sities that separates the Ising transition for naa < 1 
from the Lifshitz transition encountered at weak cou- 
pling, nas 3> 1. The corresponding effective theory has 
not been identified so far. 
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